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Abstract. We study directed last passage percolation on the planar square lattice whose 
weights have general distributions, or cquivalcntly, queues in series with general service distri- 
butions. Each row of the last passage model has its own randomly chosen weight distribution. 
We investigate the limiting time constant close to the boundary of the quadrant. Close to the 
j/-axis, where the number of random distributions averaged over stays large, the limiting time 
constant takes the same universal form as in the homogeneous model. But close to the x-axis 
we sec the effect of the tail of the distribution of the random environment. 



1. Introduction 

This paper studies the limit shapes of some last passage percolation models in random en- 
vironments, specifically, the corner growth model and two Bernoulli models with different rules 
for admissible paths. 

We introduce the corner growth model through its queueing interpretation. Consider service 
stations in series, labeled 0,1,2, ... ,£, each with unbounded waiting room and first-in first-out 
(FIFO) service discipline. Initially customers 0, 1, 2, . . . , k are queued up at server 0. At time 
t = customer begins service with server 0. Each customer moves through the system of 
servers in order, joining the queue at server j + 1 as soon as service with server j is complete. 
After customer i departs server j, server j starts serving customer i + 1 immediately if i + 1 
has been waiting in the queue, or then waits for customer i + 1 to arrive from station j — 1. 
Customers stay ordered throughout the process. Let X(i,j) be the service time that customer i 
needs at station j, and T(k,i) the time when customer k completes service with server £. 

Asymptotics for T(k, t) as k and I get large have been investigated a great deal in the past 
two decades. A seminal paper by Glynn- Whitt [6] studied the case of i.i.d. {X(i, j)}. They took 
advantage of the connection with directed last-passage percolation given by the identity 

(1.1) T(M) = max V X(i,j). 

7T & J 

The maximum is taken over non-decreasing nearest-neighbor lattice paths ir C I 2 , from (0, 0) 
to (k,£) that are of the form n = {(0,0) = (x Q ,y ), (xi,yi), . . . , (x k+ e,yk+e) = (k,£)} where 
(xi,yi) — (xi_i,7/i_i) = (1,0) or (0,1). A quick inductive proof of (jl.l|) together with earlier 
references to this observation can be found in [6] (see Prop. 2.1). This particular last passage 
model is also known as the corner growth model. 
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Next we add a random environment to this queueing model. The environment is a sequence 
{Fj : j <E Z+} of probability distributions, generated by a probability measure- valued ergodic 
or i.i.d. process with distribution P. Given the sequence {Fj}, we assume that the variables 
{X{i,j)} are independent and X(i,j) has distribution Fj. In the queueing picture this means 
that the service times {X(i,j) : i 6 Z + } at service station j have common distribution Fj, and 
at the outset the distributions {Fj : j £ Z + } themselves are chosen randomly according to some 
given law P. Obviously the labels "customer" and "server" are interchangeable because we can 
switch around the roles of the indices i and j. 

The asymptotic regime we consider for T{k,£) is the hydrodynamic one where k and £ are 
both of order n and n is taken to oo. Under some moment assumptions standard subadditive 
considerations and approximations imply the existence of the deterministic limit 

*(x,2/) = lim n~ 1 T([nx\, \ ny\) for all (x, y) € M+. 

n— >oo 

Only in the case where the distributions Fj are exponential or geometric has it been possible 
to describe explicitly the limit VP. This is the case of • /M/l queues in series, which in terms 
of interacting particle systems is the same as studying either the totally asymmetric simple 
exclusion process or the zero-range process with constant jump rate. For rate 1 i.i.d. exponential 
{X(i,j)} the limit ^(x,y) = {^fx + y/y) 2 was first derived by Rost [17] in a seminal paper on 
hydrodynamic limits of asymmetric exclusion processes. The random environment model with 
exponential Fj 's was studied in pQ [T2J [20] ■ 

Let us now set aside the queueing motivation and consider the last-passage model on the 
first quadrant Z?_ of the planar integer lattice, defined by the nondecreasing lattice paths and 
the random weights {X(i,j)}. For the queueing application it is natural to assume the weights 
nonncgative, but in the general last-passage situation there is no reason to restrict to nonnegativc 
weights. 

The ideal limit shape result would have some degree of universality, that is, apply to a broad 
class of distributions. Such results have been obtained only close to the boundary: in [14] Martin 
showed that in the i.i.d. case, under suitable moment hypotheses and as a \ 0, 

(1.2) $(l,«) = fj + 2(T\/5+o(^), 

where [i and a 2 are the common mean and variance of the weights X(i,j). Also, the o(y / a) 
term in the statement means that lim Q ^ a^ 1 / 2 [^(1, oi) — fi — 2ay/a] = 0. In the i.i.d. case ^ 
is symmetric so the same holds for &(a, 1). 

Our goal is to find the form Martin's result takes in the random environment setting. VP is no 
longer necessarily symmetric since the distribution of the array {X(i,j)} is not invariant under 
transposition. So we must ask the question separately for ^(l,a) and ^(a, 1). 

It turns out that for ^f(a,l), where the number of rows stays large relative to the number of 
columns, the fluctuations of the environment average out to the degree that our result in Theorem 
12.21 below is essentially identical to Martin's result in the homogeneous environment. We still 
have ^(a, 1) = /i + 2a^/a + o(y/a) as a \ 0, where now /x is the mean as before but a 2 is the 
average of the "quenched" variance. That is, if we let /j,q = J x dF (x) and a 2 = J(x~^ ) 2 dF (x) 
denote the mean and variance of the random distribution Fo, and E expectation under P, then 
H = EOuo) and a 2 = E(a 2 ). 

The case ^(l,a) does not possess a clean result such as the one above. Even though we are 
studying the deterministic limit obtained after n has been taken to infinity, we sec an effect 
from the tail of the distribution of the quenched mean fio. We illustrate this with the case of 
exponential {Fj}. Now the number na of distributions Fj is small compared to the number n of 
weights X(i, j) in each row, hence the fluctuations among the i^'s become prominent. The effect 
comes in two forms: first, the leading term is no longer the averaged mean /x but the maximal 
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mean. Second, if large values among the row means m are rare, the order of the a-dependent 
correction is smaller than the \J~ol seen above and this order of magnitude depends on the tail 
of the distribution of [Iq. As an exponent characterizing this tail changes, we can see a phase 
transition of sorts in the power of a, with a logarithmic correction at the transition point. 

For general distributions we derive bounds on "5(1, a) that indicate that in the case of finitely 
many distributions the correction is of order y/a. 

As auxiliary results we need bounds on the limits for last-passage models with Bernoulli 
weights under a random environment. However, with Bernoulli weights the standard corner 
growth model is not one of the explicitly solvable cases. The model with Bernoulli weights does 
become solvable when the path geometry is altered suitably. The model we take up is the one 
where the paths are weakly increasing in one coordinate but strictly in the other. There are 
two cases, depending on which coordinate is required to increase strictly. If we require the x- 
coordinatc to increase strictly then an admissible path {(xq, yo), (xi,yi), . . . , (x m ,y m )} satisfies 

(1.3) x l+ i -Xi = l and y < yi < • • • < y m . 

The other case interchanges x and y. These cases have to be addressed separately because the 
random environment attached to rows makes the model asymmetric. The sum of these two 
last-passage values gives a bound for the case where neither coordinate is required to increase 
strictly in each step. 

We derive the exact limit constants for Bernoulli models with both types of strict /weak paths. 
For one of them this has been done before by Gravner, Tracy and Widom [SJ. Their proof utilizes 
the fact that the distribution of T(k,£) is a symmetric function of the environment (at least for 
the particular Bernoulli case they study). Our proof is completely different. It is based on the 
idea in [T§] where the limit for the homogeneous case was derived: the last-passage model is 
coupled with a particle system whose invariant distributions can be written down explicitly, and 
then through some convex analysis the speed of a tagged particle yields the explicit limit of the 
last-passage model. 

Further remarks on the literature. The present paper does not address questions of fluctua- 
tions, but let us mention some highlights from the literature. For the last-passage model with 
i.i.d. exponential or geometric weights, the distributional limit with fluctuations of order n 1 / 3 
and limit given by the Tracy- Widom GUE distribution was proved by Johansson [TO]. As for 
the shape, universality has been achieved only close to the boundary, by Baik-Suidan [5] and 
Bodineau-Martin [3]. 

Fluctuations of the Bernoulli model with strict /weak paths and homogeneous weights were 
derived first in |llj and later also in [7]. For the model in a random environment fluctuation 
limits appear in [9j |8] . 

On the lattice we can imagine three types of nondecreasing paths: (i) weak-weak: both 
coordinates required to increase weakly, the type used in (jl.lj) : (ii) strict- weak: one coordinate 
increases strictly, as above in (|1.3[) : and (iii) strict-strict: both coordinates increase strictly so 
an admissible path {(x , y ), (xi,yi), . . . , (x m ,y m )} satisfies x < ■ ■ ■ < x m and y < ■ ■ ■ < y m - 
As mentioned, with Bernoulli weights the strict-weak case is solvable but the weak-weak case 
appears harder. The third case, strict-strict, is also solvable with Bernoulli weights. The shape 
was derived in [18] and recent work on this model appears in [5]. 

Organization of the paper. The main results on the shape close to the boundary are in Section 
[2] and the results for Bernoulli models in Section [3j Section [4] sketches the proof of the existence 
of the limiting shape, a result we basically take for granted. The main proofs follow: in Section 
[5] for Theorem 12.21 on ^(a, 1), in Section [6] for Theorem 12.31 on \&(l,a), and in Section [7] for 
Theorem 12.41 for the exponential model. 
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Some frequently used notation. We write 

ess sup / = inf{s eM:P(/> s) = 0} 

p 

for the essential supremum of a function / under a measure P. Z+ = {0,1,2,...}, N = 
{1, 2,3,...}, and R+ = [0, oo). 1(A) is the indicator function of event A. 

2. Main results 

First a precise definition of the last-passage model in a random environment. Let P be a 
stationary, ergodic probability measure on the space Mi(R) z+ of sequences of Borel probability 
distributions on R. E denotes expectation under P. For some of the main results P will be 
assumed to be an i.i.d. product measure. A realization of the distribution- valued process under 
P is denoted by {i r j}jez + - This is the environment. Given {Fj}, the weights {X(z) : z G Z+} are 
independent real- valued random variables with marginal distributions X(i,j) ~ Fj for (i,j) G 
Z 2 , . Let (CI, F, P) be the probability space on which all variables {Fj,X(i,j)} are defined, and 
denote expectation under P by E. 

A (weakly) nondecrcasing path is a sequence of points zq = (xo,yo), Zi = (x%, y±), . . . , z m — 
(x m , y m ) in Z 2 , that satisfy x < xi < ■ ■ ■ < x m , y Q < y x < ■ ■ ■ < y m , and \x l+ i - x l \ + \y i+ i - 
yi\ = 1. For Zi,Z2 G 1? + with Z\ < Z2 (coordinatewise ordering), let H(zi,Z2) be the set of 
nondecreasing paths from z\ to z%. Whether the endpoints Z\ and Z2 are included in the path 
makes no difference to the limit results below. The last-passage time T(z\,Z2) from Z\ to Z2 is 
defined by 

T(z\,zn) = max > X(z). 
»sn( Zl , 22 )^ 

When zi=0 abbreviate U(z) = 11(0, z) and T(z) = T(0, z). 
Put these three assumptions on the model: 

(2.1) E|A(z)|<co, 

r°° r 1 1/2 

(2.2) J |l-E(F (a?))| < oo, 

and 

/■OO 

(2.3) / csssup(l — F (x)) < oo. 

Jo p 

We begin with this by now standard result that defines our object of study, namely the function 
The proof is briefly commented on in Sectional 

Proposition 2.1. Assume F is ergodic and satisfies (|2.1[) , (|2.2p and (|2.3p . Then for all (x,y) G 
(0, oo) 2 the last passage time constant 

(2.4) V(x,y)= lim -T( [nxj , [ny\ ) 

exists as a limit both P -almost surely and in L X (P). Furthermore, ^(x,y) is a homogeneous, 
concave and continuous function on (0, oo) 2 . 

Assumption (|2.2j) is also used for the constant distribution case, see (2.5) in [TJ]. Some 
further control along the lines of assumption (|2.3[) is required for our case. For example, suppose 
1 — Fj(x) = e~^ x for random £j G (0,oo). Then (|2.3|) holds iff ess inf P (£ ) > 0. If the distribution 
of £o is not bounded away from zero, n~ 1 T(n,n) — > oo because we can simply collect all the 
weights from the row with minimal £j among {£o, ■ • ■ ,£n}- However, assumption (|2.2[) can be 
satisfied without bounding £o away from zero. 
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Now we turn to the main results of the paper on the form of the limit shape at the boundary. 
As explained in the introduction, for ^(a, 1) we find a universal form as a \ 0. In addition to 
the earlier assumptions, we need similar control of the left tail of the distributions: 

(2.5) / (E[F (x)]) 1/2 dx < oo 

J — OO 

and 
(2.6) 



,o 

/ ess sup Fq(x) dx < oo. 

J-oo P 



Let us point out that (|2.2p and (|2.5p together guarantee E|X(z)| 2 < oo. Let fij = fJ,(Fj) and 
(Tj = a 2 (Fj) denote the mean and variance of distribution Fj. These are random variables under 
P with expectations ji = E(/zo) and a 2 = E(oq). 

Theorem 2.2. Assume the process {Fj} is i.i.d. under P, and satisfies tail assumptions (|2.2|l . 
(j2~3l) . ([23]) and (|2l)jl . T/ien, as a|0, *(a,l) = + 2o\/a" + o(Va~). 

Assumptions (|2.2p and (|2.5p are direct counterparts of what was used for Theorem 2.4 in 
[14] . Assumptions (|2.3p and (|2 .6[) are additional assumptions needed for handling the random 
environment. These assumptions are used to control estimates that come from bounding limits 
of Bernoulli models. 



We turn to the case ^(l,a). The results will be qualitatively different from Theorem | 
The leading term will be the essential supremum of the mean instead of the averaged mean and 
we will see different orders for the first a-dependent correction term. 

First a general result for which we restrict ourselves to the case of finitely many distributions, 
but we can relax the i.i.d. assumption of the random distributions. 

Theorem 2.3. Assume the process {Fj} of probability distributions is stationary, ergodic, and 
has a state space of finitely many distributions Hi , . . . , Hi, each of which satisfies Martin's [14] 
hypothesis 

poo i>0 

(2.7) / (I - H^x)) 1 / 2 dx + H e (x)^ 2 dx < oo. 

JQ J -co 

Let n* = max? ^{Hi) be the maximal mean of the Hi 's. Then there exist constants < c\ < 
Ci < oo such that, as a j 0, 

(2.8) n* +c 1 y/a + o(y / a ) < *(l,a) < n* + c 2 Va + o(y/a ). 



We would expect c\ = c 2 but our proof does not give it. 

Finally, we consider the case ^(1, a) for the exponential model where some (partially) explicit 
calculation is possible. Here we see how the tail of the random mean (1q creates different orders 
of magnitude for the a-dependent correction term. Let {£j}jez + be an i.i.d. sequence of random 
variables that satisfy < c < £j with common distribution m. To distinguish the exponential 
model from the general one we write Gj(x) = 1 — e~^ jX for the distribution function of the 
exponential distribution with parameter £j, and for the limiting time constant. We assume 
c is the exact lower bound: m[c, c + e) > for each e > 0. Then the essential supremum of the 
random mean is = c . 

An implicit description of the limit shape was derived in [20] by way of studying an exclusion 
process with random jump rates attached to particles. We recall the result here. One explicit 
shape is needed for the proof of Theorem 12.21 also, so this result will serve there too. 
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Define first a critical value u* = f, ^ G (0, oo]. For < u < u* define a — a(u) 

implicitly by 

1 -m(df). 



a(u) is strictly increasing, strictly concave, continuously difFerentiable and one-to-one from < 
u <u* onto < a < c. We let a(u) = c for u > u* if u* < oo. Then define g : R + -> M+ by 

(2.9) = sup{-j/u + a(u)}, 2/ > 0. 

u>0 

The function g is monotone decreasing, continuous, and g{y) = for y > a'(0+) = 1/fJ.c- It is 
the level curve of the time constant. The equations connecting the two are g(y) = inf{a; > : 
®G(x,y) > 1} and 

(2.10) * a (x, y) = inf{f > : tg{y/t) > x}. 

Qualitative properties of the limit shape depend on the tail of the distribution m at c+, 
and transitions occur where the integrals J, — c)~ 2 m{d£) and Jj c ^{S, — c) _1 m(d£) blow 

up. (For details see [20].) These same regimes appear in our results below. For the case 
Jjc 00) — c ) 2 m (^0 — 00 wc ma kc a precise assumption about the tail of the distribution of 
the random rate: 

m[c, £) 



(2.11) 3 f e [-1, 11, K > such that lim ■ 

The value v = — 1 means that the bottom rate c has probability to{c} = k > 0. Values i/ < — 1 
are of course not possible. 

Theorem 2.4. For the model with exponential distributions with i.i.d. random rates the limit 
has these asymptotics close to the x-axis. 
Case 1: J, ^(C — c)~ 2 m(d£) < 00. Then there exists ctQ > such that 

(2.12) *g(1,£*) =c- 1 +a [ -^—m{di) forae[0,a }. 

J[c,oo) 5 — c 

Case 2: (|2.1ip /Wcfc so £/iat, in particular Jj c ^(C — c ) 2 rn{d£) = 00. Then as a \, 0, 

(2.13) i/z/ G (0, 1] tVien * G (l,a) = c _1 +a / m(d£) + o(a) ; 



'[c,oo) £ ~ c 

(2.14) if v = then ^g(1j a ) — c 1 — K ct loga + o(a log a) ; 

(2.15) if v & [—1, 0) i/ien $ G (l,a) = c _1 +Ba^+o(a^). 



In statement (|2.15p above i? = -B(c, K, v) is a constant whose explicit definition is in equation 
(|7.7[) in the proof section below. The extreme case v = — 1 is the one that matches up with 
Theorem [231 

For some heuristic understanding of Theorem 12.41 wc turn to the queueing interpretation 
discussed in the Introduction. Quantity n^!(l,a) represents the time when customer n departs 
from server \na\ (rigorously speaking in the n — > 00 limit), when initially all customers are 
queued up at server 0. When u* < 00 (Case 1 and subcase v > from Case 2) the results of [T] 
suggest that, at some distance but not too far from the first queue, the system should converge 
to an equilibrium where the queue length at server j (whose service rate is £j ) is geometric with 
mean c/ (£j — c) and the departure process from each queue has rate c. A customer arriving at a 
queue with this geometric number of customers present spends on average time 1 / (£, — c) at that 
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queue. In Case 1 this picture is precise enough so that equation (|2.12j) can be naively understood 
in these terms: a single customer travels through [na\ servers in time na f, s m(d£). After 
this it takes another n/c time to see n customers go through server IjiaJ- Together these terms 
make up the right-hand side of (|2.12[) . This argument requires a high density of customers 
because once the customer density drops below u* the system chooses an equilibrium with flow 
rate below c. (Again, for details we refer to [T].) 

This point can be detected in the proofs for Case 2 where we find a parameter clq that in some 
sense represents a flow rate and replaces c in (|2.12l) (sec cqn. (17.21) ). The formulas in Case 2 are 
then obtained by estimating c — ag. 



3. Bernoulli models with strict-weak paths in a random environment 

This section looks at last-passage models with Bernoulli-distributed weights. The environment 
is now an i.i.d. sequence {pj}j^i + of numbers pj £ [0,1], with distribution P. Given {pj}, the 
weights {X(i,j)} are independent with marginal distributions P(X(i,j) = 1) = pj = 1 — 
P(X(i,j) = 0). We consider two last-passage times that differ by the type of admissible path: 

for z±, z 2 € Zi. 

(3.1) T^(z 1 ,z 2 ) = max VlW and T*(z 1 ,z 2 ) = max Vl(z), 

In terms of coordinates denote the endpoints by Zk = (afc,6fc), k = 1,2. Then admissible 
paths 7T e n^(zi,z 2 ) are of the form tt = {(ai,y ), (ai + l,Vi), (ai + 2,j/ 2 ), • • ■ , (a 2 ,y a2 ~a 1 )} 
with b\ < yo < yi < ■ ■ ■ < y a2 ~ai < b 2 , while paths tt € H^(zi,z 2 ) are of the form n — 
{(xo, bi), (xi,bi + 1), . . . , (xi, 2 _6 15 b 2 )} with ai < xq < x% < ■ ■ ■ < x b2 - bl < a 2 . Thus paths in 
n_j.(zi,Z2) increase strictly in the rr-direction while those in Tl^(zi,z 2 ) increase strictly in the 
y-dircction. The last-passage times T^(z\, z 2 ) and T-j- (zi, z 2 ) record the maximal weights of such 
paths in the lattice rectangle ([ai,a 2 ] x [61,62]) H r L\. The following figures illustrate the two 
types of admissible paths when we take z\ — (0, 0) and z 2 = (5, 5): 



3 + 3 




(a) (b) 



Figure 1. Admissible paths in Il^(zi 1 z 2 ) and H^(zi, z 2 ) 
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As before we simplify notation with T_>.(0, z) = T^(z). The almost sure limits are denoted 

by 

(3.2) = lim -T^([nx\,[ny\) and ^ t (x,y) = lim -T t ([nx\ 7 [ny\) 

for € (0,oo) 2 . The existence of the limits needs no further comment. 

The next theorem gives the explicit limits. (|3.3[) is the same as in [§J Thm. 1]. Inside the 

E( ■ • ■ ) expectations below p is the random Bernoulli probability. Let b = ess sup p denote the 

p 

maximal probability. 

Theorem 3.1. The limits in (|3 . 2[) are as follows for x, y € (0, oo). 

( bx + y(l-b)E(^), x /y>E(f^) 
(3-3) = V&ij^r) V, E(t^) < x/y < E(f^) 

[x, 0<x/y<E{ J ^) 
with zq £ (b, 1) uniquely defined by the equation 

p(l -p) 



x/y = E 



(3.4) 



■(^o-p) 2 - 



V, x/y>E(i^) 
with zq € (0, oo) uniquely defined by the equation 

p(l-p) 



x/y = E 



{zo+p) 2 



Our second result gives simplified bounds that are useful for the proof of the main result 
Theorem 12.21 Let p = E(p) be the mean of the environment. ty(x,y) is the limiting time 
constant with weakly increasing paths defined in Proposition 

Theorem 3.2. The following three inequalities hold for the Bernoulli model: 

(3.5) 



®^(x,y) < bx + 2y/p(l - b) 



xy, 



(3.6) 

and 

(3.7) 



*t( x ) y) <pv + 2i/p(i -p)xy 



^(x, y) <py + 4y / p(l - p)xy + bx. 



(|3.7|) follows from (|3.5|l and (|3.6|) because *f?(x,y) < ^^,{x,y) + ^^(x,y). Another loose 
estimate we will use later following (|3.7[) is 

(3.8) *(.t, y)<py + iy/p(l - p)xy + bx<(y + 4y/xy)y/p + bx. 

We prove the formulas and inequalities first for and then for W-j-. For some parts of the 
proofs it is convenient to assume b < 1. Results for the case 6=1 follow by taking a limit. 
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Proof of p.3j) and (|3.5|) . We adapt the proof from [19] to the random environment situation 
and sketch the main points. 

Consider now the environment {pj} fixed, but the weights X(i,j) random. For integers 
< s < t and a, k define an inverse to the last passage time as 

T((a,s),k,t) = mm{l G Z+ : T_>((a + 1, s + 1), (a + l,t)) > k}. 

Note that T((a, s), 0, t) = but L((a, s),k,t) > for fc > 0. Knowing the limits of the variables 
r is the same as knowing *&_>.. By the homogeneity of it is enough to find h(x) = fy^(x, 1). 
By the homogeneity and superadditivity of 'J'-*., h is concave and nondecreasing. Let g be the 
inverse function of h on K + . Then g is convex and nondecreasing, and 

tg{x/t) = lim -r((0,0), \nx\, \nt\). 

n— foo Ji 

To find these functions we construct an exclusion-type process z(t) = {zk(t) : k G Z} of 
labeled, ordered particles z/c(i) < Zk+i(t) that jump leftward on the lattice Z, in discrete 
time t G Z+. Given an initial configuration {^(0)} that satisfies z,_i(0) < Zi(0) — 1 and 
liminfi-^-oo |i| _1 Zi(0) > —1/6, the evolution is defined by 

(3.9) z k (t) = inf {zi(0)+T(( Zi (0),0),k-i,t)}, k G Z, t G N. 

It can be checked that z(£) is a well-defined Markov process, in particular that Zk(t) > — oo 
almost surely. 

Define the process {rji(t)} of interparticlc distances by r).i(t) — z i+ i(t) — z.i(t) for i e Z and 
t G Z_|_. By Prop. 1 in [19j process {r]i(t)} has a family of i.i.d. geometric invariant distributions 
indexed by the mean u G [1, b~ ) and defined by 

(3.10) P(r H =n) = u- 1 (l-u- 1 ) n - 1 , n G N. 

Let Xfe(i) = Zfc(t — 1) — Zfc(t) > be the absolute size of the jump of the fcth particle from time 
t — 1 to i, and let g< = 1 — p t . From (6.5) in [19], in the stationary process 



(3.11) P(x k (t)=x) = 



(1 - up t )q t 1 a; = 

p t (l - up t )q^ 1 (u-l) x (uqt)- x x = 1,2,3, 



We track the motion of particle z (t) in a stationary situation. The initial state is defined by 
setting z (0) — and by letting {%(())} be i.i.d. with common distribution (|3.10[) . With k = 0, 
divide by t in (|3.9[) and take t — > oo. Apply laws of large numbers inside the braces in (|3.9I) . 
with some simple estimation to pass the limit through the infimum, to find the average speed of 
the tagged particle: 

(3.12) — lim -^o(^) = sup{ux — g(x)} = f(u). 

The last equality defines the speed / as / = g + , the monotone conjugate of g. It is natural to 
set f(u) = for u G [0, 1), fib' 1 ) = /((6" 1 )-), and f{u) = oo for u > b' 1 . By [M Thm. 12.4] 

(3.13) g(x) ~ sivp{xu — g + (u)} = sup {xu — f{u)}. 

u>0 l<u<l/b 

Since zq{£) is a sum of jumps xo(k) with distribution (|3.11|) we have the second moment bound 
sup 4gN E[(t _1 ^o(0) 2 ] < °°; an( i consequently the limit in (|3.12p holds also in expectation. From 
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this 



/(«) 



lim F,[t~ 1 z (t)] = lim eIV 1 Vx (fc) 



(3.14) 



fe=i 



E^x(u-l) x (uq)- x p(l-up)(l-p)- 1 = E 



E[x o (0)] 
pu(w — 1) 



up 



We find 5 (x) from (l3~T4")l and (I3TT31 : 



a: > 6 2 Eii^£l - 1 

— (b-p) 2 

1-p (b-p)- 
< a; < Ej^ 



where u Q 6 (1, b^ 1 ) is uniquely defined by the equation x + 1 = E(l — p)(l — u p)~ 2 . From this 
we find the inverse function h(x) = g^ 1 (x) and then ^_>.(x, y) = yh{x/y). We omit these details 
and consider (|3.3p proved. 

To prove f|3. 5|) we return to the duality (|3 . 1 3[) and write 



(3.16) 



g(x) > sup {xu - f(u)} for f(u) 

1<m<1/6 



u(u — 1) 
1 -ub 1 



f'(u) = x is solved by u* = b 1 (^1 — 

When ie > we have u* g [1, -|), and then 



6a: +p 



1 2 

> XU* - f{u) = —(yJ{l- b)p- y/bx+p) . 



Consequently 



g -\x) < {Vte+vtr^r-z 



bx-p + 2^(1 - b)px. 



When x < -^-r, the supremum in (|3.16[) is attained at u = 1, and in this case 



g^ix) <x<bx + 2y/(l - b)px. 
The bound (|3.5[) now follows from \I'_ ) .(x, y) = yg~ 1 (x/y). 



□ 



Proof of (|3.4[) and (|3.6[) . The scheme is the same, so we omit some more details. The inverse of 
the last-passage time is now defined 

T((a,s),k,t) = min{Z g Z+ : T t ((o, s + 1), (a + l,t)) > k}. 

Vertical distance t — s allows for at most t—s marked points, so the above quantity must be set 
equal to oo for k > t—s. The particle process {z(t) : t g Z + } is defined by the same formula ()3.9[) 
as before but it is qualitatively different. The particles still jump to the left, but the ordering 
rule is now Zk{t) < Zk+i(i) so particles are allowed to sit on top of each other. Well-defmedness 
of the dynamics needs no further restrictions on admissible particle configurations because the 
minimum in (|3.9p only considers i g {fc — t, . . . , k} so it is well-defined for all initial configurations 
{z t (0) such that z t (0) < z i+ i(0). 

The following can be checked. Under a fixed environment {pj}, the gap process {rji(t) = 
Zi+i(t) — Zi(t) : i g Z} has i.i.d. geometric invariant distributions P{r]k = n) = (i^jj)(iq3jj)™j 
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n G Z+, indexed by the mean u £ R+. In this stationary situation the successive jumps Xk(t) 
Zkit — 1) — Zfc(t) of a tagged particle have distribution 



P(x k (t)=y) 



l+up t 



y = 



l_u_\y Pt 



> 1. 



From here the analysis proceeds the same way as for the other model. The speed function is 
defined by 



f(u) = - lim E[n- 1 z (n)] = E[x o (0)] = u(u + 1)E 



P 



1 + up 



and then convex analysis takes over. We omit the remaining details of the proof of (|3.4p 
To prove (|3.6p . note that 



r \ r ti w ^> / pu(u + l)\ 

(x) = swp{xu — }{u)\ > sup<^ xu — > 

u >0 »>n ^ 1 + up J 



u>0 ' 



pu(u + 1) 

up 



I < x < p. 



We used Jensen's inequality and concavity of p > 1+ p . From this 



'p-px + 2^p(l-p)x 0<x<±-fi 



x > 



and (13.61) follows. 



□ 



4. Proof of Proposition 12. II 

We comment briefly on the proof of Proposition [2TT1 Further details can be found in [T3]. The 
flow of arguments is standard. First one takes an integer point (x, y) € and applies Liggett's 
version of the subadditive crgodic theorem to the process Z m ^ n = —T((mx,my), (nx,ny)), < 
m < n, to prove that ^f(x,y) exists and is finite. Then rational (x,y) and real (x,y) are 
handled by approximations. Along the way regularity properties of ^ are established and used: 
supcradditivity, homogeneity, concavity and continuity. 

All this works easily for the Bernoulli case because last-passage times are uniformly bounded 
in terms of path length. Consequently we can assume that Proposition 12.11 has been proved for 
the Bernoulli case. For the general case we check that for integer points (x, y) € the moment 
hypotheses of the subadditive ergodic theorem [4] p. 358] follow from our assumptions (|2.1[) . 
(ET21 . and (12311 : 

EZ o + 1 <E|T((0,0),(x, 2 /))|<E £ \X(i,j)\ = (x + + 1)E|X(0, 0)| < oo. 

0<i<x,0<j<y 
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Next: 



— EZo n > E max 

n 11 Tr£n(nx,ny) 



1 f 

V X(z)+ = -E max V / I(X(z) > u) du 

z£k 1 ' y 'z£lT J{ > 

1 f°° 1 f 00 

> E / max I(X(z) > u) du = / E max I(X(z) > u) du 

n Jn w&Tl{nx,ny) — ' n Jr. ireUtnx.ny) — ' 

>-/ sup-E max y^I(X(z)>u)du = - V Ber[1 _ F ( u)] (x,y) du 

Jo n n weU{nx,ny)^^ J 

>-(y + iy/xy) / v/l -EF (u)du-x / (l - essinfF (u))< 
Jo ./0 p 



7(yl) is the indicator function of event A. ^ Ber[L-F(u)]{ x \ y) is the limiting time constant for 
the Bernoulli model where the weights have distributions P(X(i,j) = 1) = 1 — Fj(u) = 1 — 
P(X(i,j) = 0). On the last line above we used the Bernoulli estimate (|3.8|) . By assumptions 
(|2.2j) and (|2.3I) . EZo,„ > nj for a constant 7 > —00. These estimates justify the application of 
the subadditive crgodic theorem. We omit the remaining details and consider Proposition 12.11 
proved. 



5. Proof of Theorem 

For the first lemma, let {Fj} and {Gj} be ergodic sequences of distributions defined on a 
common probability space under probability measure P. In a later step of the proof we need 
to assume {Fj} i.i.d. Assume that both processes {Fj} and {Gj} satisfy the assumptions made 
in Theorem 12.21 With some abuse of notation we label the time constants, means, and even 
random weights associated to the processes {Fj} and {Gj} with subscripts F and G. So for 
example fip = E(J xdFo(x)). The symbolic subscripts F and G should not be confused with the 
random distributions Fj and Gj assigned to the rows of the lattice. We write <5 Ber([G(x)-F(x)] + ) 
for the limit of a Bernoulli model with weight distributions P(X(i,j) = 1) = (Gj (x) — Fj(x))+ 
I — P(X(i,j) = 0) where x is a fixed parameter. An analogous convention will be used for other 
Bernoulli models along the way. 

Lemma 5.1. Assume {Fj} and {Gj} satisfy ([2~2]) . (|2"T3j) . (|2~5|) and ([2~6l> . Then for a > 0, 



(5.1) 



|* F (a, 1) - *Gf(a, 1) - (p F - Hg)\ 

f+00 i/ 2 n+00 

< 8^ / (E|G (x) -F (x)\J dx + a esssnp\F (x) - G (x)\dx. 



Proof. The right-hand side of (|5.1j) is finite under the assumptions on {Fj} and {Gj}. Couple 
the Fj and Gj distributed weights in a standard way. Let {u(z) : z = £ be i.i.d. 

Uniform(0, 1) random variables. Set Xp(z) = F^ 1 (u(z)), where Fj~ (u) = sup{a; : Fj(x) < u}, 
and similarly Xg{z) = Gj (u(z)). Write E for expectation over the entire probability space of 
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distributions and weights. 

-* G (a,l) 

= lim — E max ~S^X F (z) — lim — E max X G (z) 

moon 7ren(l cm.J .ri) Z — ' n-s-oo rt iren( I anj ,n) z — ' 



n-toofl 7ren(LanJ,n) ^ 



< lim -E max V(X F (z) - X G (z)) 

n— >OC n ttCT\(\ ntm I ml ^ ' V ' 

/■+OC 

= lim -E max V / I l(X G (z) <x< X F (z)) - l(X F (z) < x < X G (z)) \dx 

nosn 7r6n(laril,n) ^ /_„ I V ' V 

< lim — E / max S^{l(X G (z) < x < X F (z)) - l(X F (z) < x < X G (z)) Xdx. 

n^oo n J_ ca n£ll(lan\,n) ^ I V ' V 7 J 

We check that Fubini allows us to interchange the integral and the expectation. Since F and G 
are interchangeable it is enough to consider the first indicator function from above. Let a be an 
integer > a. 

r +00 



-E max l(X G (z) < x < X F (z))dx 

_ x n 7ren(Lanj,n)^ v w; 

/+00 ^ 
sup -E max l(X G (z) < x < X F (z)) dx 
„ n 7ren(on,n) J 

*Ber([G(x)-F(x)] + )(a, 1) 

< / ' (e|G (x)- F (x)\ + 4^a~(E\G {x) ~ F (x)\) 1/2 + aesssiip\G a {x) ~ F {x)\J dx 



—00 
+00 



00 
+00 



< 00 



by estimate (|3.7p and the finiteness of the right-hand side of (|5.ip . Continue from the limit 
above by applying Fubini. Then take the limit inside the dx-integral by dominated convergence, 
justified by the n-uniformity in the bound above. Finally apply again the Bernoulli estimate 

y F (a,l)-V G (a, 1) 

: lim / ^E max V(/(X G (z) < x < X F (z)) - l(X F (z) < x < X G (z)) X dx 
lim — i E max V l(X G (z) < x < X F (z)) 

+ E ^ X( 1 - / (*H*)<*<*<K*)))-X;i}<fc 

/+00 
{*Ber([G(x)- F (x)]+)(a, 1) + ^ Ber(l-[F(x)-G(x)] + ) («, 1) ~ (1 + «)}dx 
-00 

< / + (e(G (x)-F (.t)) + + 1-E(F (x)-G (x)) + 



+ 4 v / q t ( ^E(G (x) - F (x)) + + yjE(F (x) - G (x))_ 
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+ a ( ess sup [Go (x) — F a (x)] + + 1 — ess inf [F (x) — Gq(x)] + ) — (1 + a) >dx 
\ p p J J 

/+oo /-+00 
\/E\F a (x) — Gq(x) \ dx + a / esssup|Go(x) — F (x)\ dx. 
-oo J-oo P 



Interchanging F and G gives the bound from the other direction and concludes the proof. □ 

For a while we make two convenient assumptions: that the weights are uniformly bounded, 
so for a constant M < oo, 

(5.2) F{F (-M) = and F (M) = 1} = 1, 

and that variances are uniformly bounded away from zero, so for a constant < cq < oo, 

(5.3) P{<7 2 (F ) > c } = 1. 
Note that then 

(5.4) c < cr 2 (F ) < M 2 P-a.s 

and the conditions assumed for Theorem 12.21 arc trivially satisfied by the uniform boundedness. 

Henceforth r = r(a) denotes a positive integer-valued function such that r(a) / oo as a \ 0. 
Tile the lattice with 1 x r blocks B r (x, y) = {(x, ry + k) : k = 0, 1, r — 1} for (x, y) G Z 2 , . A 
coarse-grained last-passage model is defined by adding up the weights in each block: 

x r ( Z )= Yl x ( v )- 

veB r {z) 

The distribution of the new weight X r (i, j) on row j S Z + of the rescalcd lattice is the convolution 
F r j = F r j * F r j + i * ■ ■ ■ * F r j +r -\. 

We repeat Lemma 4.4 from P3] with a sketch of the argument. 

Lemma 5.2. Let typ(x,y) and ^F r (x,y) be the last passage time functions obtained by using 
Fj and F r _j as the distributions on the jth row, respectively. If r — > oo and ry/a as a 1 0, 
then 

lim — j=\^ F (oi, 1) '3/ p (ar, 1)1 =0. 

aio y/a ' r r ' 

Proof. Given a path it £ Il(m, nr — 1), consider all the blocks that it intersects; this gives a path 
7f £ II(m, n — 1) in the rescaled lattice such that I {U z( z% B r (z)) /\tt\ < mr. Then by (|5.2[) 



max y j X(z) — max N ^(z) <mrM. 



Take m = [anr\ , divide through by nr, and the conclusion follows. □ 
Let and V r ^ v be the mean and variance of F rtV : 

r— 1 r— 1 

i=0 j=0 

Let <& r; y be the distribution function of the normal A/"(/i r ,y, V r , y ) distribution, and $, vy the 
distribution function of Af(rfip, V r . y ). The difference between $> r ,y and $ rjJ) is that the latter has 
a non-random mean. We shall also find it convenient to use {Xj} as a sequence of independent 
variables with (random) distributions Xj ~ Fj. For the next lemma we need to assume {Fj} an 
i.i.d. sequence under P. 

As in [13] , a key step in the proof is the replacement of the rescaled weights with Gaussian 
weights, which is undertaken in the next lemma. 
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Lemma 5.3. Assume {Fj} i.i.d. under F. If r — > oo and r^/a — s- as a I 0, then 



(5.5) 



lim - 

aio r-Ja ' 



|*p (or, 1) - 4-* far, 1)1 = 0. 



Proof. According to Theorem 5.17 of |15j . independent mean random variables Xi, X 2 , X3, . . . 
satisfy the estimate 



p{B-^jr X t < x] *(*) I < + | X |)-3 J , 

i— 1 B r 



x e k, 



where B r — X)i=i Var(X-j), $ is the standard normal distribution function, and A is a constant 
that is independent of the distribution functions of X\, X 2 , ■ ■ . , X r . Then. 



(5.6) 



\F r , v (x) ~ * r ,,(*)| < A ^^!^ ^ + ' |3 (1 + V-^\X - » r J) 

\ r 



-3 



where the second inequality used the assumptions < M) = 1 and of > c§ > 0. 

Armed with (|5.6p we now estimate the right-hand side of (|5 . 1 [) for the processes {F r ^y\ y ^%, 
and {^ > r,y}yez + and with a replaced by ar. 

For the first term on the right in (|5.1[) , note this Schwarz trick: for a probability density / on 
R and a function H > 0, 



y V~Hdx = J f^y/f-iHdx < (jf^Hdx 



1/2 



For the calculation below take <5 > and /(a;) = ci(l + |x — r/ip\ 1+s ) 1 for the right constant 
ci = ci((5). Factors that depend on M and 6 are subsumed in a constant C. Then 



oar f (E|F r , (x) - $ n o(z)|) 



1/2 



f/.r 



< CaVW / < E 



(l + M-V-^-Mr.ol)" 



1/2 



(/;)' 



KCaVWl® (1 + 
I J —00 

by a change of variables x = /U r + yM,Jr 



^ii+^-^W] 1/2 



< Ca 



1 / 2 r 1 /2| E |^ _ r(UF |l+5 + r d+5)/2| 1/2 < C . a l/2 r (3+5)/4^ 

In the last step we used E|// ri o — r(ip\ 1+s < Cr^ 1+S '^ 2 which follows because fx r ,o — r/iF is a sum 
of r bounded mean zero i.i.d. random variables. 
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For the second term on the right in (|5.1D . 

esssup|F r (x) — <f> r o(x)\ dx < Car 1 / 2 / esssupllH -=^- I dx 

-co V J-oo p V M V r J 

„ 1/2 f f~™ ( -rM-x\~ 3 , r rM , /"+ 00 / i-rMV 3 , 

< CW 1/2 <M 1 + =- cfe + / dx + / 1 + =r- 

U-oc V M^F J J-rM JrM V ^ 

< Car 3 ' 2 . 

To summarize, with these estimates and (|5.1[) we have 

J-|* F >r, 1) - **>r, 1)| < JL( a VV 3 + 5 )/ 4 + ^3/2^ 



If (5 is fixed small enough, assumptions r — > 00 and r^/a — ► make this vanish as a — >• 0. □ 

The next lemma makes a further approximation that puts us in the situation where all sites 
have normal variables with the same mean. 

Lemma 5.4. Let *f/$ r and be defined as before, and again r\fa — > as a — > 0. Then 

]im—=\^^(ar,l)-^ (cw,l)| = 0. 
aio ry/a r 

Proof. For z = G l\, let have distribution $ rj so that = X^ r \z) - fj, rJ + 

r\iF has distribution & r ,j- Now estimate: 

*j (or, 1)= lim - max y^J? (r) (z) 

r n-i-oo n TrenXLcmrJ ,n) Z — ' 

< lim — max 'S^X^(z) + lim — max /x r + r/x F ) 

n-i-oo n 7rGn(LQnrJ .n) Z — / n-i-oo 77, 7rGll( [cmrj ,n) Z — ' 

1 " 1 

< flj*.(ar, 1) + lim — N (— fj, r j + rfip) + lim — 2Mr ■ I anr\ 

r n— >oo 71 ^ — ' v ' n— ±00 77, 

= *$ r (ar, 1) + 2Mar 2 . 

The opposite bound (ar, 1) > \E'$ i .(ar, 1) — 2Mar 2 comes similarly, and the lemma follows. 

□ 

Let us separate the mean by letting <j> ry denote the N(0,^2lZo a ry+i) distribution function. 
Since the last-passage functions of the normal distributions satisfy (ar, 1) = r/j F (l + ar) + 
fy-^( r )(ar, 1), we can summarize the effect of the last three lemmas as follows. 

Lemma 5.5. Assume {Fj} i.i.d. under V, and assume r = r(a) satisfies r — > 00 and r^/a — > 
as a I 0. Under assumptions (|5.2|) and (|5.3j> 

(5.7) lim -^=\9 F (a, 1) - - -*s« K, 1)1 = 0. 



In order to deduce a limit from (|5.7[) we utilize the explicitly computable case of exponential 
distributions from |20j . We need to match up the random variances of the exponentials with 
the variances a 2 of the sequence {Fj}. Thus, given the i.i.d. sequence of quenched variances 
a 2 = a 2 (Fj) that we have worked with up to now under condition (|5.4j) . let £j = 1/oj and 
Gj{x) = 1 — e~^ x the rate £j exponential distribution. Then {£j}jez + is an i.i.d. sequence of 



LAST-PASSAGE SHAPE 



17 



bounded random variables < c < £j < b with distribution m. We can assume c is the exact 
lower bound: m[c, c + e) > for each e > 0. Gj has mean and variance ^{Gj) = £~ and 
a a (G J )=^ 2 = 4 

Assumptions (|2.2p and (|2.3p are easily checked, and so the last-passage function is well- 
defined. We would like to apply Lemma l5.5l to this exponential model, but obviously assumption 
(|5.2[) is not satisfied. To get around this difficulty we do the following approximation which 
leaves the quenched means and variances intact. We learned this trick from |14j . 

Let Yj denote a Gj-distributed random variable. For a fixed r > 0, let 

mj = EiY^Yj > t) and w j =E(Y?\Y j > r). 

The quantities 

K~ T ) 2 and u _ w i - 7-2 _ 

— - — 7J ana u> — T 

[rrij — T) z + wj — rrij — r 

satisfy the equations 

(f — Sj)r + sjUj = rrij and (1 — Sj)r 2 + SjU 2 = Wj. 
Then < sj < 1, Uj > t and Wj > t 2 . Define distribution functions 

!Gj(a;) < a; < r 

1- Si [l-Gj(r)] r<x<Uj 
1 X > Uj. 

Yj — Gj satisfies Bij = E'Y, and i?!^ 2 = EY 2 . Moreover, for any fixed r > 0, 

_ E(Y 2 \Y 3 >r)-r 2 _2 2 

so the distributions {Gj} are all supported on the nonrandom bounded interval [0,2/c + r\. 
Consequently Lemma 15.51 applies to G. We can draw the same conclusion for G once we have 
the next estimate: 

Lemma 5.6. Given e > 0, we can select t large enough and define Gj as in (|5.8p so that 

1 

alO y/a 



lim— |*g(Q!j 1) - *g(«, X )l < £ - 



Proof. This comes from an application of Lemma l5.f I Gj = Gj on (— oo,r) and 1 — Gj < f — Gj 
on all of R. The integrals on the right-hand side of (|5.f [) are finite and can be made arbitrarily 
small by choosing r large. □ 

Currently we have shown that 
(5.9) lim^|* G (a,f)-Ea --% (or,l)| =0. 

It remains to perform an explicit calculation on ^Q(a, 1). As before, utilize the notation /i G — 
E^ 1 and a 2 G = E^ 2 . 

Lemma 5.7. For random exponential distributions with rates bounded away from zero, 

* G (a, 1) = fJ>G — 2a G y/a~ + 0(a). 
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Proof. Recall the definition of the limit shape ^g{<^,1) from (|2.10j) . From ()2.9[) one can read 
that tg(l/t) is nondecreasing in t. Thus by (|2.10p ^g(oc, 1) = t = t(a) such that tg(l/t) = a. 

Next we argue that when a is close enough to 0, g(l/t) = —uo/t + a(uo) for some < uq < u* 
with a' (ito) = 1/t. Since a(0) = and a(u*— ) = c, strict concavity gives for < u < u* 



...(#)} 1 = a'(u*-) <o'(«) = { / f m(dO} 1 

<o'(o+) = {/ r'mCdo)' 1 



e,oo) J Mg 

On the other hand, < ^gC^i 1) — MG < C\/a + Ca where the second inequality comes from 
comparing {Gj} in (|5.1j) with identically zero weights. Thus when a is small enough, 1/t is in 
the range of a'. Consequently there exists uq £ (0,u*) such that a'(uo) = or equivalcntly, 

(5.10) f m(de) = *. 

From the choice of t, a = tg{\/t) = t[—ua/t + a(i*o)) = — + ta(uo) and so 

a ur\ a f 1 

(5.11) t G a, 1 = t = — — + = — — + / p^rdm 0- 

a(u J a(u ) a(u ) J[ c ,oo) ? - a(«oJ 

Combining (|5.10[) and (|5.11[) gives 

(5.12) a = a(uo) 2 / y- — --^m(d£). 



From this 



/ -H m (dO < a. 

J[c,oo) C 



Hence we have < a(ito) < y/a/ac- 

When a and hence a(uo) is small, f|5.12[) and the last bound on a(u ) yield 



^ "T"^2 - G G 



f 1 1 

/ [?7 — „r„ fia ~ tsM^) 
i[ c ,oo) (.? ^ a{u )) 4 



(5-13) = / ^K) aM 2 w 

3) ru-a(uo))^ 



Consequently 



77 Ks mW) = 0(V5). 

o) ?(c- a(u )) 2 



a G = — = 0{y/a). 



a(u ) + o~ G 
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Now we put all the above together to prove the lemma. 



* G (a, 1) - pt G - 2a G ^/a 
a f 1 



a(u ) 

a 
a(u ) 



i{d£,) — /j,g — 2cr Gy /a 



+ -^a(u ) + O(a(u ) 2 ) m(d£) - /XG - 2a G ^/a 



) £ - a(u ) 
1 1 

( + ( 



a(u ) 
0(a) as a I 0. 



og ) + cr G a{u ) (a G 



a(u ) 



□ 



Combining Lemma l5~7l and (|5.9[l gives 

lim I (ar, 1) - 2cr G v/al = 0. 
aio y/a r r 

Substitute this back into (|5.7[) and recall that a F = a G . The conclusion we get is 

(5.14) lim-^|* F (a, 1) - \i F -2a F ^/a\ = 0. 

We have proved Theorem 12.21 under assumptions (|5.2I) and (|5.3I) . We now lift (|5.3[) . For 
£ > 0, let {W^(z)} be i.i.d weights with distribution H defined by P(W(z) = ±e) = 1/2. Let 
Fj = Fj * H be the distribution of the weight X(i, j) = X(i, j) + W(i, j). Let ^>h and \E'^ be the 
time constants of the last-passage models with weights {^(z)} and {X (z)}, respectively. The 
Bernoulli bound p.7[) gives the estimate ^n{x 7 y) < Ae^/xy. The corresponding last-passage 
times satisfy 

Tp[z) - T H {z) < T F (z) < Tp(z) + T H {z) 
where Tjj(z) uses the weights —W(z). In the limit 

(5.15) ® F (a, 1) - 4eV« < 1) < 1) + 4e\/a- 

Since a 2 (Fj) = a 2 (Fj) + e 2 while /Us = fj, F , and e > can be arbitrarily small, this estimate 
suffices for limit (|5.14|) . 

As the last item of the proof of Theorem [22] we remove the uniform boundedness assumption 
(|5.2|) . Suppose {Fj} satisfy the conditions required for Theorem 12.21 but there is no common 
bounded support. For a fixed M > define the truncated distributions 

x > M 

—M <x<M 
x < -M. 

Let hm, &m and ^ FM {x,y) be quantities associated to {-Fj,m}- 
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From (|5.1j) and the conditions assumed in Theorem 12.2 



< 



^E|F (a;) - F 0>M (x)\j dx + ^fa j esssup\F (x) - F QtM (x)\dx 

r 

(E|i^ (x)| 



v 1/2 

!|1-F (x)|) dx 



-M 



-t-oo 



esssup|Fo(x)| dx + / esssup|l — F (x)\ dx 
'mp 



< e. 



The last inequality comes from choosing M large enough, and is valid for all a < 1. Since 
E(i?J 2 (0,0)) < oo, dominated convergence gives om — >• a and so we can pick M so that 
|c — crjvf | < e. Now 

-=|tf F (a!,:L)-M-Wa| ^ -^|*F M (Q!,l)-/iM-2a- Af v / a|+2e. 

Since e is arbitrary and limit (|5.14[) holds for {Fj^m}, we get the conclusion for the sequence 
{Fj}. This concludes the proof of Theorem 12.21 



6. Proof of Theorem 



Proof of Theorem \2.3[ The lower bound in (|2.8[) can be proved by applying Martin's result (|1.2[) 
to the homogeneous problem where a maximal path is constructed by using only those rows j 
where Fj = if.;*, the distribution with the maximal mean (j,* = u(Hi*). This is fairly straight- 
forward and we leave the details to the reader. 

To prove the upper bound in (|2.8p . we start by increasing all the weights X(z) by moving 
their means to n* . Then we subtract the common mean \i* from the weights, so that for the 
proof we can assume that all distributions H\, ... , Hl have mean zero. 

Create the following coupling. Independently of the process {Fj}, let {Xg(z) : 1 < I < L, z £ 
r i?j r } be a collection of independent weights such that Xg(z) has distribution Hi. Then define 
the weights used for computing ^(1, a) by 



X(z) =Y,X t {z)I {Fj=Hl} for z = G 



Begin with this elementary bound: 



= lim -E 

n— yoo Ti 



(6.1) 



max y X(z) 

Trenfn.lcmJ) ^ 



< lim — E max Xp(z)Isp.—ij,\ 



The next lemma contains a convexity argument that will remove the indicators from the last 
passage values above. 

Lemma 6.1. Let T> be a sub-a-field on a probability space (tt,J-,P), D an event in T>, and £ 
and r\ two integrable random variables. Assume that En = 0, rj is independent of T>, and £ and 
r\ are independent conditionally on T>. Then E[£ V (j]Id) ] < E[£ V n]. 
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Proof. By Jensen's inequality, for any fixed 

x V E(r) | V) < E{x\J T}\V). 
Since r\ is independent of T> and mean zero, 

iV0< E{x\lr\\V). 

Integrate this against the conditional distribution P(£ £ dx\T>) of £, given 2?, and use the 
conditional independence of £ and 77: 

V0|X>) < E^VtjIV). 

Next integrate this over the event D c : 

E[l D c- £V{r)I D )] =E[lDo-tV0] <E[Idc- £Vti]. 
The corresponding integral over the event D needs no argument. □ 

Fix a lattice point zq = (io,jo) for the moment. We split the maximum in (|6.ip according to 
whether the path tt goes through zq or not, and in case it goes we also separate the weight at zq: 

n m ¥ ,,52 X t(*)hFi=n t } = (B + X e (z )I {F He} ) =B+(A-B)V (X e (z )I {Fja=He} ) 
where 

A = ma,xy^X e (z)I< F H} and B = max X t (z)Is F i= H e }- 

z£k\{z } 

Now apply Lemma T6. II with £ = A — B, r\ = X^(zo); an d -D = {Fj = Hg}. Given Fj , A — B 
does not look at X^zq), so the independence assumed in Lemma |6. II is satisfied. The outcome 
from that lemma is the inequality 



E 



max Vl((z)I(,,H ( ) <E[AV (B + X e (z ))]. 

7ren(n, [an] ) '-^ 



This is tantamount to replacing the weight Xn(z§)I{ Fj =H e } a t z q with Xg(zo). 

We can repeat this at all lattice points zq in (|6.1|) . In the end we have an upper bound in 
terms of homogeneous last-passage values, to which we can apply Martin's result (|1.2p : 

l ^ L 

*(l,a)<V lim -e[ max XAz)] = V (1, a) 



7rGn(n, [cm] ) 



fcl J/ zStt £=1 

L 

= 2Va53<r(fli) +o(V5). 

This completes the proof of Theorem 12.31 □ 

7. Proofs for the exponential model 
Proof of Theorem \2.J\ Equation (|2.10l) gives 

(7.1) * G (1, a) = M{t > : tg(a/t) > 1} = t(a) = t. 

That the infimum is achieved can be seen from (|2.9I) . 

Under Case 1 the critical value u* = f, < c(£ — c) _1 m(d£) < 00, and also 
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By the concavity of a and ()2.9j) . for < y < a'(u*—) we have g(y) = —yu* + c. Consequently 
for small enough a 

1 



1 = tg(a/t) = —ac / m(d£) + ct 

J[c,oo) ? — c 

and equation (|2.12[) follows. 

In Case 2 a'(0+) > a'(u*— ) = and hence for small enough a > there exists a unique 
ito G (0,u*) such that a'(ito) = a/t. Set cio = a(«o) € (0,c). As a \ 0, both u /* u* and 
a c. We have the equations 

/" £ * 

a'(uo) -1 = / j- r^m(d£) = -, 1 = tg(a/t) = -au + ta , 

j[ C ,oo) (£-aor a 

(7.2) *g(1,«) = * = - + — = - + a / -J— m(d£) 

a a a J[ c ,oo) 4 _ a o 

and 

1 /" 1 

(7.3) — = a 77 T^MO- 

Assuming (|2.1ip . start with i/ € (— 1, 0) U (0, 1). For a small enough e > there are constants 
< K\ < k.2 such that 

(7.4) Kl (£_ c )"+i < m [c,0 <K 2 (C-cr +1 br(e[c,c + e] 
and as e \ we can take Ki, ki k. First we estimate c — ao- Fix e > 0. 

for a quantity Ci(er) = 0(e~ 2 ). The first term above can be bounded above and below by (|7.4I) . 
and we develop both bounds together for Ki, i = 1,2, as 

fc+e ft _ n \v+l 



1 u+1 



r +£ [(e-tto)-(c- a o)]' 

2nia — — ^ d£ + Ci(e) 



fv 4- 1 \ /" c+e 

= 2^5^ ; )(-i) k (c-a ) k (z-a Q y- k - 2 dt+c 1 ( E ) 

(7.5) fe=o V K 7 Jc 

= 2 ^«oE t (-l) fc (c-a )^ . '- + d(e) 



fc=0 



k J k—u+1 
2K ia 2 A u (c - ao)"- 1 - 2 Ki a 2 ,, ^ ; (« - «o) fe (c + e - ao)"-*- 1 + tt(e) 



fc / fc - z/ + 1 

fc=0 v 7 

= 2 Kia 2 y l,( c - ao )^- 1 +C 1 (e). 
Ci (e) changed of course in the last equality. In the next to last equality above we defined 

fe=0 



k I k-v+l 
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Rewrite the above development in the form 

(c - ao) 1 ^ = 2nc 2 A u a + a[2A„(«jfl$ - kc 2 ) + Ci(e)(c - ao) 1- *]. 

Now choose e = e(a) \ as a \ but slowly enough so that Ci(e)(c — ao) 1_ly — > as a \ 0. 
Then also and wc can write 

i i 

(7.6) c — ciq ~ B^a 1 -" + o{a 1 - 1 ') 

with a new constant Bq = (2kc 2 ^4^) t ^ 77 . 

Now consider the case v £ (0, 1) which also guarantees Jj c — c) _1 m(d£) < oo. From (|7.2p 
and (T776I) as a \ 



1 /" 1 

* G (l,a) = — +a / m(d£) 

a J [coo) ? - a 

= -+a [ —^—m(d£) + 0(aT J ^) + a([ — - — m(d£) - / — !— m(d£) 

c ,/[c,oo)C;- c V ./[coo) ? ~ a J[c,oo) ? ~ c 

1 /" 1 

= - + a m(d£) + o(a). 

c 7[ c ,oo) ? - c 

Next the case v e (—1, 0). The steps are similar to those above so we can afford to be sketchy. 

1 f 1 

* G (l,a) = ha/ m(d£) 



a J[c,oo) £, — a 

1 , c~a Q , (c-a ) 2 , ^[CO , ~ M 



c c 2 c 2 a J c (£ - a ) 2 

Again, using (|7.4[) and proceeding as in (|7.5p . we develop an upper and a lower bound for the 
quantity above with distinct constants /Cj, i = 1,2. After bounding m[c, £) above and below with 
— c) !y+1 in the integral, write (£ — c)' y+1 = ((£ — a ) — (c — a )) l/+1 and expand in power 
series. 

1 i i f c+e (£ - c) 1/+1 
- + %" 2 q 1 -" +o(ai-") + aK, / ^j- J —- T d£ + aC 1 (e) 



c J c (C-ao) 2 

1 _o -J— , -J_n / s„ (v + l\ (-l) fe 

= - + B c a 1 -- + o(a ) + aK,(c - a ) 2^ 



/c I k — v 

k=0 x 7 



fe=0 

1 

c 



v — k \c— ao + e 
+ BaT^" + o(qi^) + A„ )2 a(«i - «)(c - ao)" + aCi(e) 



In the last equality the next to last term with the X)fc°=o sum was subsumed in the aCi(e) term. 
Then wc introduced new constants 

(7.7) ^,2 = X)( I )Vt^ and B = B c- 2 + K B»A^ 2 . 

k=0 ^ / V 

As before, by letting e = e(a) \ slowly enough as a \ we can extract ^0(1,0;) = c _1 + 
Bq 1 ^ + o(a T ^ 77 ) from the above bounds. 

It remains to treat the cases v = —1,0,1 where integration of the type done in (|7.5[) is 
elementary. We omit the details. □ 
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